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Abstract 

We investigated the viscoelastic properties of colloid-polymer mixtures at intermediate colloid 
volume fraction and varying polymer concentrations, thereby tuning the attractive interactions. 
Within the examined range of polymer concentrations, the samples ranged from fluids to gels. 
Already in the liquid phase the viscoelastic properties significantly changed when approaching the 
gelation boundary, indicating the formation of clusters and transient networks. This is supported 
by an increasing correlation length of the density fluctuations, observed by static light scattering 
and microscopy. At the same time, the correlation function determined by dynamic light scattering 
completely decays, indicating the absence of dynamical arrest. Upon increasing the polymer con- 
centration beyond the gelation boundary, the rheological properties changed qualitatively again, 
now they are consistent with the formation of colloidal gels. Our experimental results, namely the 
location of the gelation boundary as well as the elastic (storage) and viscous (loss) moduli, are 
compared to different theoretical models. These include consideration of the escape time as well 
as predictions for the viscoelastic moduli based on scaling relations and Mode Coupling Theories 
(MOT). 

PACS numbers: 62.20.-x,62.10.+s,64.70.pv 
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I. INTRODUCTION 



The mechanical properties of sohds and hquids are very different. A sohd responds 
elastically to a small deformation, while a liquid flows. Complex fluids, such as colloidal 
suspensions, polymers or surfactant solutions, have mechanical properties between those of 
elastic solids and viscous liquids, they are viscoelastic. Furthermore, they can be significantly 
perturbed by even modest mechanical forces. This causes a wealth of fascinating effects 
l|, but also provides a challenge to fundamental and applied research to understand their 
behavior under deformation and flow, i.e. their rheology. A detailed knowledge of their 
properties is crucial for many applications; complex fluids are extensively used in industrial 
products and processes Q, 

Among complex fluids, colloidal suspensions are frequently used as models of atomic 
systems whose interparticle interactions can be tuned 4|. These model systems allow the 
investigation of various fundamental phenomena, such as the equilibrium thermodynamics 
of gas, liquid and crystal phases and also the non-equilibrium behavior of gels and glasses. 
The latter are disordered solids which are dynamically arrested and long-lived. Although 
they can be formed at any colloid volume fraction , most studies have focused on either 
very large or small volume fractions. 

At low colloid volume fractions, the interparticle attraction induces the formation of 
clusters 0, y, 3 which may interconnect to create a space-spanning network U 
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A connection between the gelation boundary and the spinodal line has been proposed 
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J. 

m 

16l . Il7l | with spinodal decomposition driving cluster formation and gelation [l^ . The 



elasticity of gels is related to the connectivity of the network and the size of the clusters, 
i.e. to the heterogeneous structure of the network 10, Is], 19, 20]. 

In contrast, at large volume fractions amorphous solids are alr eady formed in the absence 



of attraction due to crowding; repulsive hard sphere glasses 
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22l ]. Increasing attraction 



strength initially causes melting of the glass, before for even larger attractions again amor- 



phous solids, attractive glasses, form [23|, l2J, l25| . Mode-coupling theory (MCT) has shown 



that the formation of repulsive and attractive glasses is caused by dynamical arrest due to 



caging and bonding, respectively 
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281 ]. These two arrest mechanisms seem to 



dominate also the mechanical response to deformations. The elastic properties are deter- 
mined by the confinement of particles and can be rationalized in terms of the ratio between 
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the energy and the volume characterizing the structural length of the system, which at high 
volume fractions coincides with the particle size HQ. I. addition, «ea,a.d scat- 



tering experiments on repulsive 
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32| and attractive glasses 
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3J| indicated one (in the 



former) or two (in the latter) step yielding related to cage and bond breaking. 

In the region of intermediate colloid volume fractions the origin of the fluid-solid transition 
(even in the quiescent state) is still under debate and the rheological properties of the 
amorphous solid are hardly understood. It has been suggested by experiments and MCT 
calculations 
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36l | that the route to gelation is comparable to the one at higher volume 



Tactions, i.e. dynamical arrest. This is, however, in contradiction to recent simulations 
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37( 1 which predict an important role of spinodal decomposition also at intermediate 



volume fractions. Moreover, it is not clear whether elasticity is caused by similar mechanisms 
as in attractive glasses or by the connectivity of the network, as suggested by the structural 

and confocal microscopy experiments 



heterogeneity observed in scattering 
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41|. 



Finally, the evolution of the elastic response inside the gel has been rarely studied. It has 
been recently suggested that the micromechanic response of gels with different interparticle 
attraction is strongly correlated with the 'clusterlike' or 'stringlike' nature of the gel structure 
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Here we are interested in this intriguing region of intermediate colloid volume fraction. In 
particular, we investigate the transition from liquid-like to solid-like behavior upon increasing 
interparticle attraction. We use a mixture of nearly hard-sphere colloidal particles with 



non-adsorbing linear polymer 



0,0,0, 



44 



451 ]. The polymer induces a depletion attraction 



between the particles whose range and strength can be tuned by the polymer size and 
concentration, respectively. We investigate the static and dynamic properties of samples with 
an increasing polymer concentration, i.e. increasing strength of the attractive interaction. 
Upon increasing polymer concentration, the system evolves from an equilibrium liquid to a 



non-equilibrium, dynamically arrested gel 



.e syst 



471 ]. We use static and dynamic light 



scattering and microscopy to investigate this liquid-solid transition and determine structural 
parameters, such as the characteristic correlation length of density fluctuations, and dynamic 
properties, such as the collective dynamics. The samples are also subjected to shear and 
their mechanical response is determined by rheology. This combination allows us to relate 
the static and dynamic properties in the quiescent samples to the behavior under small 
deformation, i.e. the linear rheological properties. We particularly focus on the evolution of 
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the viscoelastic moduli as a function of increasing attraction, i.e. on moving from liquid-like 
to solid-like samples and deeper into the gel region. Viscoelastic measurements approaching 
gelation are compared to different models: estimates of the escape time and scaling relations 
as well as MCT predictions for the time and polymer concentration dependence of the shear 



moduli. MCT has on 
shear 
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y recently been extended to describe dynamically arrested states under 
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II. MATERIALS AND METHODS 



A. Samples 

We investigated mixtures of Polymethylmethacrylate (PMMA) colloids and linear 
Polystyrene (PS, from Polymer Laboratories) dispersed in cis-decalin at a temperature 
T = 23°C. The hydrodynamic radius of the PMMA particles, R = 137 nm, was determined 
by dynamic light scattering in the very dilute regime. The polydispersity of the colloids 
was not directly measured, but suppression of crystallization in quiescent and sheared col- 
loidal dispersions indicates a polydispersity of about 12%. The radius of gyration of the PS 
(molecular weight M^, = 132.9 kg/mol) in cis-decalin, = 10.8 nm, was estimated based 



on 



571] and the polydispersity was cited as M^/M„ = 1.01. In dilute solution, this implies a 



polymer-colloid size ratio ^ = Vg/R = 0.079 ± 0.013. The effective polymer-colloid size ratio 
takes into account the concentration dependence of the polymer size and the mesh size 
in the semidilute reg ime. It has been calculated according to the Generalized Free Volume 



Theory (GFVT) |58|,l59|. 



The colloid stock solutions with volume fraction = 0.6 were prepared by redispersing 
spun-down sediments, whose volume fraction was estimated to be </> = 0.67 when taking 
polydispersity into account 60|]. Polymer stock solutions were prepared by adding cis- 
decalin to dry polymer. Polymer concentrations Cp (mass/volume) were calculated from 
the weighed masses of the two components and their densities. Colloid-polymer mixtures 
were obtained by mixing appropriate amounts of colloid and polymer stock solutions. After 
mixing, samples were vigorously shaken using a vortex shaker, then homogenized over 3 
days in a rotating wheel mixer. The compositions of the samples are summarized in table [T] 
with the first column giving the nominal polymer concentration which is used to refer to 
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TABLE I: List of samples, is the colloid volume fraction, Cp/c* and "^/c* are the polymer 
concentrations in the total and free volume, respectively, in units of the overlap concentration, 
is the effective polymer-colloid size ratio. 



samples in the following. The polymer overlap concentration c* has been estimated by 



c* = 3Mw/4:7iNAr^- Values of and Cp refer to the 



occupied by colloids as estimated by GFVT 



58 



;otal volume and c!l''^® to the volume not 
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61|. 



Samples were mixed and the light scattering or microscopy measurements perfomed 
within an hour to avoid effects due to aging (Sec. Ill Pp . 



B. Light Scattering 

Due to the difference in refractive index n between PMMA and cis-decalin (tipmma = 1.49, 
rzdec = 1-48) the samples are turbid. Multiple scattering was suppressed and single scattered 



light recorded using a 3D dynamic light scattering instrument (LS Instruments) 62|, |63 |. 
From the cross-correlation functions we extracted the dynamic structure factors /(Q,t) 
with the delay time r, the modulus of the scattering vector Q = (4ndcc7r/A) sin(6'/2), the 
scattering angle 6 and the wavelength A = 633 nm (HeNe laser from JDS Uniphase). 

In static light scattering experiments, the Q-dependence of the time-averaged intensity 
{I{Q)) was calculated from the time-averaged intensities recorded by the two detectors, 
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{Ia{Q)) and {Ib{Q)), and the intercept Pab of the cross-correlation function: 

{im = ^{i^'\Q)){ii'\Q)) = ^{ia{Q)){h{Q)mt/f3!:^) (1) 

where the superscript '(1)' refers to quantities determined in the single-scattering regime. 
To achieve ensemble averaging, the sample was rotated continuously. Rotation does not 
affect Pab, {Ia{Q)) and {Ib{Q)), but the time dependence of fiQ.r). Static structure factors 
S{Q) were obtained from {I{Q)) taking into account the particle form factor (as determined 
in the dilute regime) and the transmissions Tq and Tf, of the sample: 



where the superscript '(d)' refers to quantities determined in the dilute regime. This assumes 



that all significant contributions to the scatterin g ar e due to the colloids, as shown in 64 1 



for similar colloid-polymer mixtures (see also 65l. l66l|). 



C. Microscopy 

DIG microscopy experiments were perfomed using a Nikon Eclipse 80i upright microscope 
with a Nikon 100 x Plan Apo objective and a Canon EOS 30-D digital camera. Samples 
were loaded into a home-built cell: Two nr. 1 coverslips were glued onto a microscope slide, 
leaving a 3 to 4 mm wide channel between them. The channel was filled with the sample 
and a further nr. 1 coverslip used to cover the sample at the top and glue (UV-cure adhesive. 



Norland Optical Adhesive (NOA) 61) to seal the open ends of the channel [67 1. 



D. Rheology 

We used an ARES-HR rheometer with a force balance transducer lOFRTNl and a cone- 
plate geometry (cone angle 0.044 rad, cone diameter 25 mm) which provides a constant 
strain throughout the sample. The geometry surfaces were mechanically roughened to avoid 
wall slip. To test reliability of the geometries with roughened surfaces, we compared results 
obtained with roughened and smooth surface geometries for samples where the presence 
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of wall slip could be excluded. The agreement was found to be satisfactory. In dynamic 
measurements wall slip apparently has no effect in the linear viscoelastic regime, but dra- 
matically affects measurements at large strains in the non-linear regime when the polymer 
concentration is comparable or larger than c*. Here we only investigate the linear regime 
with strain amplitudes 0.001 < 70 < 0.02, while the results in the non-linear regime are 
discussed elsewhere 68 



In order to minimize solvent evaporation, a solvent saturation trap was used. The trap 
isolates the sample from the surrounding atmosphere by a fluid seal at the top and a per- 
manent seal at the bottom. Solvent evaporation leads to a saturated atmosphere inside the 
enclosure. 

In order to eliminate the effect of sample loading and aging, the following procedure was 
adopted: After loading, a dynamic strain sweep test was performed, i.e. the samples were 
subjected to oscillatory shear at a frequency of u = 1 rad/s and the strain amplitude 70 
was increased until the sample showed a liquid-like response; 70 = 8 was sufficient at all Cp. 
Moreover, before each test, oscillatory shear with uj = 1 rad/s and 70 = 8 was imposed on the 
samples until G' and G" reached constant, steady-state values. Subsequently, samples were 
left at rest for a waiting time before the test was started. We performed aging experiments, 
a series of Dynamic Time Sweeps at = 10 rad/s and a total duration of 50000 s without 
any rejuvenation in between. Samples below the macroscopic gelation boundary showed 
no aging effects over the whole time interval, while gels showed an initial increase of the 
elastic modulus within the first 200 s after loading, but then the moduli remained constant 



at least up to 3600 s. A detai 
will be reported elsewhere 



ed study of aging effects on the rheological properties of gels 
691 ]. Here we note that for 200 s < < 3600 s the viscoelastic 
properties did not change and reproducible results were obtained in consecutive tests. We 
have chosen t,„ = 300 s. 



7 



2.0 
1.5 
1.0 
0.5 
0.0 

0.30 0.35 0.40 0.45 0.50 

FIG. 1: Macroscopic behaviour as investigated by tube inversion for different volume fractions (j) 
and polymer concentrations Cp/Cp. (•) gels (no flow), (■) highly viscous fluids, (A) low viscosity 
fluids. 



gelation boundary ^ 



III. RESULTS AND DISCUSSION 



A. Quiescent Samples 

1. Macroscopic Behavior 

The macroscopic behavior of the samples was investigated by tube inversion as a function 
of colloid volume fraction and polymer concentration Cp/c* (Fig. [1]). Gel samples were 
identified by the absence of flow after tube inversion, which was found for Cp/c* > 0.4. 
Samples with 0.2 < Cp/c* < 0.4 showed already a relatively high viscosity, but were still 
flowing. 



2. Microscopic Structure 

Microscopic structural information was obtained by static light scattering (SLS) and 
DIG microscopy. We determined the static structure factor S{Q) at low scattering vectors 
0.35 < QR < 2.6 (Fig. [2]) where length scales corresponding to collective structures of 
the order of a few particle diameters are probed. Due to the limited Q-range we cannot 
observe the first peak of S{Q) which for a pure colloidal dispersion with (p = 0.4 is expected. 
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FIG. 2: Static structure factor S{Q) measured by static light scattering for different polymer con- 
centrations Cp/c* (as indicated in the legends) for samples below (a) and above (b) the macroscopic 
gelation boundary. 



depending on polydispersity, in the range 3 < QR < 4 [6^ 

Below the macroscopic gelation boundary {cp/c* < 0.4, Fig. [2J\), S{QR<1) monotoni- 
cally increases with increasing Cp/c*. For Cp/c* < 0.25 a finite value of S((5^0) could be 
extrapolated which is consistent with the clustering of particles due to attractive deple- 
tion interactions, as has already been observed for silica-PS mixtures at the same colloid 
volume fraction 65|]. These clusters are not necessarily equilibrium clusters . For 

larger polymer concentrations, Cp/c* = 0.32 and 0.4, S{QR<1) increases steeply, which in- 
dicates an increasing amplitude of the density fluctuations. Crossing the gelation boundary 
{cp/c* ^ 0.4), S{QR<1) drops dramatically pointing at the suppression of large density fluc- 
tuations. Then S{QR<1) increases again inside the gel region, i.e. for Cp/c* > 0.4 (Fig. [2|3). 
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FIG. 3: Upper plot: Static structure factor S{Q) at QR = 0.35 (red circles) and QR = 1.4 (blue 
squares) as a function of polymer concentration Cp/Cp. Lines are fits by a power-law dependence 
for Cp/c* < 0.4, and a linear dependence for Cp/c* > 0.4. Inset: Correlation length ^ obtained by 
fitting an Ornstein-Zernike scaling to the low Q part of S{Q). Lower pictures: DIG microscopy 
images of samples with polymer concentrations < Cp/c* < 2 as indicated. 



The Cp-dependence of S{Q) at two distinct Q values {QR = 0.35 and 0.7) is summarized in 
Figure [31 The strong increase of S{Q) upon approaching the gelation boundary, i.e. Cp/c* < 
0.4, can be described by a power-law dependence, S{Q) ~ {cp/c*)°' [g^, with an exponent 
= 4.6 ± 0.3 for QR = 0.35 and ao.7 = 2.6 ± 0.1 for QR = 0.7. After the sharp drop at 
the gelation boundary, S{Q) increases roughly linearly with increasing Cp/c* inside the gel 
region. 

Fitting an Ornstein-Zernike scaling, S{Q) ~ l/[<5^ + (1/C)^] to those S{QR<1) which 
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increase at low Q, a characteristic correlation length ( can be extracted (Fig. [3l inset). We 
find that C/2-R increases from approximately 1 to 6 upon increasing Cp/c* from 0.32 to 0.4 
and then, inside the gel region, drops again to approximately 1 with a slight increase with 
increasing Cp/c*. 

Due to the limited Q-range accessible in our light scattering experiments and, as a con- 
sequence, the large uncertainty in the value of (, we complemented our light scattering 
experiments by DIG microscopy (Fig. [3]). At Cp/c* = the sample appears homogeneous 
refiecting its fiuid structure. Increasing Cp/c* toward the gel boundary, some graininess due 
to large scale structures is visible with the length scale and amplitude increasing strongly at 
the gelation boundary (cp/c* = 0.4). Within the gel phase (cp/c* > 0.7) the length scale and 
amplitude of the observed graininess decreases and subsequently saturates. DIG microscopy 
thus indicates that structural heterogeneities have a maximum around Cp/c* = 0.4 and their 
length scale well before gelation and inside the gel region appears to be comparable. This 
is consistent with our light scattering results. 

Large scale heterogeneous structures have been reported based on confocal microscopy 
experiments for gels composed of polymer-grafted silica spheres 39|] and in dense PMMA-PS 



suspensions 



71| where the maximum degree of heterogeneity has also been observed in 



correlation 
In this 



the vicinity of the gelation boundary. The sharp maximum of the structural 
length at the gelation boundary is expected for an arrested phase separation 
scenario, phase separation at the gelation boundary leads to a coarsening and cluster forma- 
tion that is interrupted by dynamical arrest when the clusters permanently bond to form a 
gel. Due to our relatively large volume fraction, it is also conceivable that clusters connect 
to form a transient percolated network, which only arrests when the bond lifetime becomes 
large enough. To confirm that phase separation ultimately causes gelation is beyond the 
scope of this report and requires a more detailed study, including an investigation of the 
time dependence of the low Q scattering after mixing. 



3. Dynamics 

Dynamic light scattering (DLS) was performed with samples below the macroscopic gela- 
tion boundary {cp/c* < 0.4). Measurements were done at different scattering vectors Q, 
all of them below the first peak of the structure factor where the dynamic structure factor 
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FIG. 4: Dynamic structure factor /(Q, r) at QR = 0.7 as a function of delay time r for samples with 
polymer concentration Cp/c* = (red), 0.1 (blue), 0.2 (green), 0.25 (turquoise) and 0.32 (violet; 
from left to right). Inset: Initial decay of /(Q,t). Lines are fits to the short-time expansion of 
f{Q,T) derived from the Smoluchowski equation j^. 



f{Q,T) reflects collective dynamics. The data obtained at QR = 0.7 are shown in figures H] 
and [5l That f{Q,T) completely decays indicates that the particle dynamics is ergodic, 
consistent with the macroscopic gelation boundary at Cp/c* > 0.32. The decay of f{Q,T) 
cannot be described by a single exponential: At short times it is exponential, at long times 
stretched exponential. With increasing Cp the particle dynamics is observed to slow down 
on both time scales (Fig. |4]). 

The initial, fast decay of f{Q, t) covers only about 3% of the total decay (Fig. HI inset). 
At short times, individual particles diffuse freely as reflected in the linear time dependence. 
This linear dependence is described by the short-time limit expression derived from the 
Smoluchowski equation j^: f{Q, r) = 1 — Ds{Q)Q'^t + O(r^), where Ds{Q) is the effective 
short-time diffusion coefficient. Departure from this free diffusion is observed at progressively 
shorter times for increasing Cp. This Cp-dependent departure is clearly visible if f{Q,T) is 
plotted as a function of the rescaled time tQ^Ds (Fig. [5l inset). The rescaling with Ds, 
whose Q-dependence is here omitted since Q is fixed, also accounts for the trivial dependence 
on the viscosity of the polymer solution rjrTjAcc {Vdcc is the solvent viscosity). 
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FIG. 5: Dynamic structure factor f{Q,T) at QR = 0.7 as a function of the rescaled delay time 
tQ^Ds with the effective short-time diffusion coefficient Ds- Samples are as in figure HI Inset: 
Initial decay of /(Q, t). 



The long-time relaxation is also slowed down (Figs.|4]and[5]). It corresponds to a collective, 
slow relaxation process related to the diffusion of particles whose movements are restricted 



Dy size polydispersity, which 
, and/or a distribution of 



by their mutual attraction. Its stretched form can be caused 
leads to a spread in the long-time self diffusion coefficients 
particle diffusivities, caused by heterogeneities in the particle density, in agreement with the 
static light scattering data (Fig. [3]). 

The short-time behavior for different Q is summarized in figure [61 Shown is the Cp- 
dependence of the normalized short-time diffusion coefficient Ds{Q)rir/DQ where Dq = 
k^T/QTfq^ccR is the free diffusion coefficient. It decreases with increasing Q and Cp, es- 
pecially in the range < Cp/c* < 2. These trends of Ds{Q)rir/ Dq = H{Q) / S{Q) result 
from a delicate balance between the structure factor S{Q) and the hydrodynamic function 
H{Q) (Fig.[6]B, C). S{Q) has been determined by static light scattering (Fig. [2]). It is almost 
constant for Cp/c* < 0.2 and then decreases toward the gelation boundary, except for the 
smallest scattering vector, QR = 0.7 (Fig. EB)- Based on Ds{Q)rir/ Dq and S{Q), H{Q) was 
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FIG. 6: (A) Normalized short-time diffusion coefficient Ds{Q)rir/ Dq, (B) structure factor S{Q) 
and (C) hydrodynamic function H{Q) as a function of polymer concentration Cp/c* for different 
scattering vectors QR = 0.7 (•), 1.35 (■), 1.9 (a) and 2.35 (♦). 



calculated (Fig. [6p). Except for the smallest Q, QR = 0.7, H{Q) decreases monotonically 
with increasing Cp. 

The values of Ds{Q)rir/ Dq for small Cp/c*, i.e. close to the pure hard-sphere case, are 
larger than 1 due to the dominant contribution of H{Q) as compared to S{Q), in agreement 
with previous data j^. Increasing attraction within Cp/c* < 0.2 hardly affects the structure, 
hence S{Q) is about constant, while it causes a decrease of H{Q) which for Cp/c* ~ 0.2 
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FIG. 7: (A) Normalized long-time diffusion coefficient DL{Q)rir/ Dq and (B) stretching exponents 
/3 as a function of polymer concentration Cp/c* for different scattering vectors QR = 0.7 (•), 1.35 
((■), 1.9 (A) and 2.35 (♦). 



becomes smaller than S{Q) resulting in a decrease of Ds{Q)f]r/ Dq to about 1. For larger 
polymer concentrations, Cp/c* > 0.2, the particle attraction starts to affect the average 
structure with a reduction of S{Q) near the peak (large QR) and an increase of S{Q) at 
small QR due to aggregation (Fig. [2]A.). Thus the trends of H{Q) and S{Q) become similar 
and Ds{Q)rir/ Dq tends to a constant value. Only for the largest polymer concentration, 
the value of Ds{Q)f]r/ Dq remains approximately constant for QR > 0.7. A similar trend 
has been observed in the liquid phase of a phase separating colloid-polymer mixture at 
comparable volume fractions 72|. 
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The normalized long-time diffusion coefficient DL{Q)rir/ Dq and stretching exponent (3 
are obtained by fitting a stretched exponential to the long-time behavior of the dynamic 
structure factor (Fig. [7]). Di{Q)rir/ Dq shows first a modest and then a sharp decrease 
with increasing Cp. This indicates that particle attraction slows down particle motions and, 
in the vicinity of the gelation boundary (cp/c* < 0.4), the particles start to be localized, 
signaling the approach of dynamical arrest. The stretching exponent (3 is always below 1 
with an approximately constant value of about 0.8 for Cp/c* < 0.25 and a sharp decrease 
at Cp/Cp = 0.32. This indicates that, upon approaching the gelation boundary, density 
fluctuations increase and broaden the distribution of effective long-time diffusion coefficients, 
consistent with the increasing correlation length ( of structural heterogeneities observed by 
static light scattering and microscopy (Fig. [3]). 



B. Samples under Shear 

Dynamic frequency sweeps (DFS) are reported in figure [S] for samples below (A) and 
above (B) the macroscopic gelation boundary, Cp/c* = 0.4 (Fig. [1]) with the frequency lo 
given in units of the instrument (top axis) and in units of the inverse diffusion time in the 
dilute limit tq = R'^/Dq ^ 4 x 10"^ s (bottom axis). 

Without polymer {cp/c* = 0), the response is characteristic for concentrated hard-sphere 



suspensions [73|: Viscous properties {G") dominate over elastic properties (G")> but the 
elasticity is still finite, i.e. G' > 0. 

For Cp/c* = 0.1 and 0.2, G" still exceeds G' over the whole frequency range. Both 
mechanical moduli show the same frequency dependence, which can be described by a power- 
law with an exponent of about 0.55. A G' which is larger than for hard-spheres arises from 
enthalpic contributions due to the interparticle attraction. The observed response shows 
interesting similarities to that measured in chemical and physical polymer gels, in par ticular 
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76(. This 



in partially cured or weakly cross-linked materials at the percolation point 
suggests that already for 0.1 < Cp/c* < 0.2 a percolated network is formed. This network 
is transient, since, in contrast to chemical gels, the lifetime of the physical bonds between 
colloids is finite with bonds dynamically forming and breaking. This dynamic, transient 
structure is consistent with an ergodic, complete relaxation of the dynamic structure factor 
/(g,r) (Fig.SD. 
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FIG. 8: Dynamic frequency sweeps with elastic (storage) modulus G' (filled symbols) and viscous 
(loss) modulus G" (open symbols) as a function of frequency a; in units of the instrument (top 
axis) and in units of the inverse diffusion time in the dilute limit tq = /Dq « 4 x 10~^ s (bottom 
scale). Samples below the macroscopic gelation boundary {cp/Cp < 0.4) are shown in the upper, 
in the gel region {cp/Cp > 0.4) in the lower plot with polymer concentrations Cp/c* = (♦), 0.1 
(T), 0.2 (♦), 0.25 (A), 0.32 (■), 0.4 (•), 0.5 (♦), 0.7 (t), 0.8 (♦), 1.0 (a), 1.5 (■), 2.0 (•). Lines 
are Mode Coupling Theory predictions for G' (solid lines) and G" (dashed lines) (Sec. [IVB|). The 
straight lines in (A) indicate the typical scaling in a newtonian liquid, G' ~ co'^ and G" ~ to. 
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The first indication of a solid-like response is found for Cp/c* = 0.25 at large frequencies, 
with the crossing point of G' and G" at ujtq = 10^^. This corresponds to structural relaxation 
times in the experimental time-window, in agreement with the fluid-like relaxation observed 
by dynamic light scattering (Fig. Hj). Due to the increasing strength of the attractions, 
the lifetime of the particle network becomes comparable to the examined time-scale (or 
frequency) and is long enough to cause solid-like behavior at short times, corresponding to 
high frequencies. 

For Cp/c* = 0.32 the frequency dependence of both, G' and G", becomes weaker and 
they cross close to the low frequency limit of the investigated frequencies. This solid-like 
response over almost all measured frequencies indicates, upon increasing Cp, an increase 
of the structural relaxation time with an increase in the lifetime of the network and the 
approach to dynamical arrest. This is consistent with the drop in the long-time diffusion 
coefficient (Fig. [7]). 

When the macroscopic gelation boundary {cp/c* = 0.4) is crossed, the frequency depen- 
dence of G' and G" is comparable at all Cp, which suggests a structural relaxation time 
consistently larger than the experimental observation time and particle dynamics which 
are arrested and hence non-ergodic samples. Thus, the percolated network lacks (measur- 
able) structural relaxation with a very long lifetime of particle bonds. Increasing Cp further 
increases the elastic response as indicated by the increase of G'. 

Interestingly, within the gel region {cp/c* > 0.4), G" presents a minimum at intermediate 



frequencies. This has alreac 
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78|, 



79 



80, 



been observed for a large variety of so-called 'soft glassy 



821 ] . Such a minimum suggests long-time (low frequency) 



structural relaxation, a relaxation, inside the gel phase. Its frequency could be associated 
with a transition from (3 to a relaxation (both outside the experimental time window) 
and thus be related to the length scale over which particles diffuse before they reach the 
transient non-ergodicity plateau between the two processes. Since the minimum stays at 
ujE? / Dq ^ 4 X 10~^ up to Cp/c* = 1 and then shifts toward lower frequencies, this indicates 
a constant a relaxation time for gels up to Cp/c* = 1 and an increasing a relaxation time 
for Cp/c* > 1. 

The Cp-dependence of the elastic modulus G' is summarized in figure [91 Three different 
regimes can be distinguished. For Cp/c* < 0.25 (region I in the figure), G' is very small 
and modestly increases with increasing Cp. (At the lowest frequency uj = 0.1 rad/s, G' 
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FIG. 9: Elastic (storage) modulus G' extracted from dynamic frequency sweeps as a function of 
polymer concentration Cp/c* at different frequencies uj = 0.1 rad/s (♦), 1 rad/s (a), 10 rad/s (■), 
100 rad/s (•). Upper inset: G' as a function of Cp/c* approaching gelation (0.25 < Cp/c* < 0.4). 
Solid lines represent power-law fits with exponents from about 6.2 (red) to about 8.3 (black). 
Dashed lines represent exponential fits. Lower inset: G' as a function of Cp/c* inside the gel phase 
{cp/Cp > 0.5). Lines represent linear fits. 

could not be detected for Cp/c* < 0.25.) For a fixed uj and increasing Cp, G' increases as 
a result of the increasing strength of attraction and the increasing entropic contribution 
from density fluctuations which are enhanced due to attraction-induced clustering. The 
density fluctuations will be averaged out at long times and will thus not contribute at low 
frequencies. Moreover, due to the low Cp the bond lifetime is short and there will be no bond 
contribution to G' at long times (low frequencies). In contrast, at higher frequencies both 
contributions are present and G' thus increases with frequency. 

Approaching the gelation boundary, 0.25 < Cp/c* < 0.4 (region II in the figure), G' shows 
a steep increase with increasing Cp (Fig. [9l upper inset), which reflects the strongly increas- 
ing number of permanent bonds and their increased strength which leads to a permanent, 
stress bearing network. The increase in G' can be described by a power-law or exponential 
dependence. The exponent of the power-law dependence increases from about 6.2 to about 
8.3 with decreasing frequency. Thus, the power law dependence tends to the exponential de- 
pendence and hence the quality of the exponential fit improves with decreasing frequency. A 
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more pronounced increase of G' with decreasing frequency is reminiscent of the discontinuos 
jump from zero to a finite shear modulus when crossing the gelation boundary as predicted 



by MCT 251]. In section HVBI we compare these data to MOT predictions, which take the 
bond energy into account, but neglect effects of heterogeneous structure and percolation. 

Within the gel region, Cp/c* > 0.5 (region III in the figure), G' increases linearly (Fig. [HI 
lower inset). The linear increase suggests that, once a gel is formed and saturation of 
permanent bonds reached, the elastic response depends on the structure of the network and 
the bond energy. A simple model which accounts for both contributions is proposed below 
(Sec UVDl) . 

The dynamic frequency sweeps and the dependence of the elastic modulus G' on polymer 
concentration Cp together with the dynamic light scattering results and macroscopic obser- 
vations indicate the existence of two transitions: First, at Cp/c* = 0.25 network formation 
and the first solid-like response with a relaxation time within the experimental time window 
is observed. Second, at Cp/c* = 0.4 gel formation and a solid-like response with a structural 
relaxation time outside the experimental time window was found. While both processes im- 
ply a network structure, the dynamics of the networks, in particular the lifetime of particle 
bonds and hence of the whole networks, seem different. 

IV. THEORETICAL MODELING 
A. Bond Lifetime 

Our results suggest that a crucial parameter is the time particles remain within the range 
of their mutual attraction, i.e. the 'bond' lifetime. This determines whether, on a given 
time scale, the network is transient or permanent. We estimate the bond lifetime as a 
function of polymer concentration Cp with a simple model. It is based on the approach by 



Kramers to describe the escape of particles from a potential well 83|]. The first passage 



time of a Brownian partic 
calculated numerically 



e within a depletion potential (Asakura-Oosawa potential) can be 
] . In order to obtain an analytical expression, we approximate the 



depletion potential by a ramp potential U (r) with the same depth Uq {Uq < 0) and width 
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2CR 

oo r <2R 

U{r) =\ Uo{l- iff) 2R<r< 2(r + l)R (3) 

r > 2(^* + 1)R 

Based on the depletion potential we estimate Uq = —IlpVo{2R), where Up is the osmotic 
pressure and Vo{r) the overlap volume of the depletion regions of two particles at distance 



58, 
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6l|, 



r. We calculated Vq according to the generalized free volume theory (GFVT) 
which accounts for the Cp-dependence of the polymer size and osmotic pressure (Sec. Ill Ap . 
The dependence of Uq on Cp/c* is shown in figure [TU] (inset). The error bars reflect the 
uncertainty in the size ratio which results from the uncertainty in the colloid and 
polymer radii. This uncertainty propagates to an uncertainty in the escape time. The 



escape time Tesc from a ramp potential is 
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84| 



D^io loo ~ d[^^ mo? 



where Dg is the short-time self-diffusion coefficient of a particle within the potential U{r). 
It is estimated based on the short-time self- diffusion coefficient of a particle in a colloidal 



dispersion having volume fraction (p = 0.4; D^^'^ ^ 0.3-Do |3[. Since the ramp potential 
overestimates the particle attraction in a depletion potential, Tesc is expected to overestimate 
the escape time for a depletion potential. 

The escape time Tesc increases rapidly with increasing Cp (Fig. [10]). At Cp/c* ~ 0.4 it 
reaches lab time scales (hours) and thus indicates permanent bonds with dynamical arrest 
and gel formation, in agreement with the macroscopic gelation boundary (Fig. [1]). For 0.1 < 
Cp/c* < 0.32, Tesc is within the experimental time window of the rheological measurements. 
In the rheology experiments we observe, within the accessible time window, a transition to 
solid-like behavior for Cp/c* = 0.25 and 0.32 (Fig. [H]). The frequency of the crossing point 
of G' and G" (Fig. [TOl filled squares) has the same order of magnitude as the calculated 
T~^. This indicates that, approaching gelation, a particle network forms whose relaxation 
time, given by the lifetime of the particle bonds, determines the structural relaxation of the 
system. For larger Cp, t^sc and the crossing point of G' and G" is beyond the time window 
accessible by rheology. For smaller polymer concentrations, < Cp/c* < 0.2, no crossing 
point was observed, although the calculated Tesc hes within the time window accessible by 
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FIG. 10: Escape time t^sc for a particle confined to a linear ramp potential U{r) as a function of 
polymer concentration Cp/c*. Tesc was calculated according to equation d] (o) and compared to the 
time at which, according to rheology experiments, G' and G" cross (■), and to the short (♦) and 
long (a) relaxation times obtained from light scattering at QR = 0.7. Dashed red lines indicate 
the time window accessible by rheology. Inset: Dependence of Uq on Cp/Cp. Line is a fit to the 
power- law dependence C/q ~ (cp/c*)'''^ 



rheology. This suggests that in these samples the structural relaxation probed by rheology is 
not related to the breaking of particle bonds and the samples are rather fluids of individual 
particles or clusters of particles than transient network structures. 

Based on the light scattering results, namely the short and long-time collective diffusion 
coefficients (Figs. EJA., [71A.), we estimate the characteristic relaxation times on a length scale 
corresponding to the range of the potential, 5\ ts = 5*'^/Df{Q) and tl = 5*^ / D^I\Q) and 
compare them to Tosc. The ratio between the short-time self diffusion coefficient, d'^'^ and 
the collective diffusion coefficient for hard spheres with (p = 0.4 as determined by dynamic 

(s) ^ 

light scattering is D's'/Ds ^ 0.2 pfl. Since attraction mainly affects the long-time decay, 
we use this ratio as an estimate for all Cp (Fig. [TOl filled diamonds). Its slight increase, 
almost invisible on the large vertical scale of the plot, is related to the Cp-dependence of 
the potential range 6*. However, this ts is much smaller than the calculated Tesc. For the 
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long-time diffusion coefficient, D^l^ /Di 0.1 for hard splieres [3j. In this case, this relation 



is only valid at small c„, since large Cp, i.e. strong attraction, slows down the 



ong-time self 



6^). This is 



diffusion due to bonding (in analogy to repulsive and attractive glasses 
consistent with the agreement observed at Cp/c* = 0, but tl increasingly underestimates the 
time needed to diffuse a distance 5*, possibly being responsible for the increasingly large 
discrepancy between tl and Tcsc when approaching the gelation boundary (Fig. [101 filled 
triangles). Moreover, we found that for Cp/c* = 0.1 and 0.2 the rheological relaxation time 
is faster than the shortest time accessible in the experiments and might thus be closer to 
Ts than Tosc. We attribute this to the fact that particles or particle clusters diffuse within 
a shorter time than the lifetime of particle bonds. This supports the above finding that 
samples with Cp/c* < 0.2 are fluids of individual particles or particle clusters rather than 
transient or permanent networks, which start to form for Cp/c* > 0.25. 

Finally, long-time structural relaxation has been observed by rheology in samples inside 
the gel region. In the gel region the bonds are so strong that they can be considered as 
essentially permanent (corresponding to huge t^sc at Cp/c* > 0.4, Fig. [10]). The long-time 
relaxation can thus not be related to bond breaking, but might be associated with particles 



85|, or with 



of different mobilities, as suggested by studies on dynamical heterogeneities 
different restructuring processes, such as rotation of particle groups, which then lead to the 
observed stress relaxation and aging. 



B. Frequency Dependence of the Moduli by Mode Coupling Theory 

Mode-coupling theory allows the shear modulus G{t) to be calculated for dense suspen- 
sions and predicts the existence of a glass transition. Within this approach, the modulus is 
approximated by calculating the overlaps of stress fluctuations with density fluctuations in 
order to capture the slow structural relaxation which occurs close to the glass transition. The 
slow relaxation of the system is thus described by the transient density correlator. While 
the full mode-couplin g eq uations possess a wavevector dependence it has been shown for 



the quiescent case j49|, |86[] that a simplified, schematic version of the theory in which the Q- 
dependence is neglected can effectively capture the essential physics. Recent developments 
generalizing the theory to the case of steady-shear have shown that a similar schematic 
model can be used to represent the full mode coupling equations under shear, the so-called 
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f1^2^ model 48|, |49|, |50|. Within this schematic model the transient density correlator (j)(t) 



obeys the equation of motion 

dt(j){t) + r (^<f){t) + ^ dt'm{t - t') {dt,(t){t') + 5 = , (5) 

where m{t) is the memory function, T is the initial decay rate and 0(0) = 1. Introduction of 
the parameter 5 provides an additional decay mechanism leading to long time relaxation of 
glassy states 86|]. The theory assumes that m(t) and the shear modulus G{t) relax on the 
same time scale as the correlator therefore a self-consistent approximation closing the 
equations of motion can be made. In the model the memory function m{t) is given by 

'^(^) = l^|^^)2 (^'l'^(^) + , (6) 

where vi and V2 are coupling vertices chosen so as to reproduce the generic behaviour of the 
full, Q-dependent theory at the glass transition and are thus not independent, but connected 



by a simple algebraic relation [48|, |49|, 150| • Typically the parameters are chosen as V2 = 2 
and Vi = f2(\/4/t>2 — 1) + e/(y^ — 1). In this way, both vi and V2 are determined by 
the separation parameter e. The value e = corresponds to the glass transition point and 
positive (negative) values of e correspond to statepoints in the glass (fluid). We note that for 
small amplitude oscillatory shear the 7 dependence in m{t) may be neglected. The modulus 
G{t) is given in the F^'^^ model by 

G{t)=vM\t)+S:S) , (7) 

where the modulus amplitude provides an additional fit parameter. The elastic (storage) 
modulus G' and viscous (loss) modulus G" are obtained by Fourier transformation 

POO 

G'{uj)+iG"{uo) = iuo dte-''''G{t)^=Q. (8) 
Jo 

In order to model the data, there are four free parameters: e, x, Wo- and F. On a double 
logarithmic plot the shape of G'{uj) and G"{uj) depends only on the distance from the glass 
transition, parameterized by e, and on the parameter x. The initial decay rate F and 
amplitude allow for horizontal and vertical translations, respectively. The additional 
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decay parameter S may then be used to fine tune G"{uj) at low frequencies. 

Tfie fundamental assumption of MCT is the description of dynamical arrest in terms of 
glassy dynamics, i.e. the slowing down of particle dynamics is caused by increasing caging of 
particles when approaching the glass transition. This limits the range of Cp, i.e. the strength 
of attraction, to which the model can be applied, namely the region approaching the gel 
transition. For low polymer concentrations, Cp/c* < 0.2, the dynamics are not adequately 
described in terms of caging. For large Cp/c*, i.e. inside the gel region, the model cannot 
reproduce the frequency dependence of G' and G" due to the different nature of the dynamics 
assumed in the model (glassy) and present in the samples. 

Predictions based on the F^^^ model for the region where it is assumed to be valid, namely 
approaching the gel transition {cp/c* = 0.2, 0.25 and 0.32), are shown in Fig. [H] (lines) . The- 
ory correctly estimates the relative magnitude of G' and G" in the high frequency limit, 
in particular for Cp/c* = 0.25, while an increasingly larger discrepancy evolves at lower 
frequencies. MCT seems to associate a Newtonian fluid response to the system after struc- 
tural relaxation, i.e. at frequencies below the crossing of G' and G", while the experimental 
response shows a less pronounced frequency dependence for both moduli. This could be 
caused by polydispersity and/or dynamical heterogeneities, which MCT does not consider 
but are indicated by the stretched exponential decay observed in dynamic light scattering 
(Sec. IIII A3|) . For Cp/c* = 0.32, G' agrees well with experiments, while the predicted G" 
shows a consistently different frequency dependence. The predicted minimum appears only 
as a shoulder (at higher frequencies than the predicted minimum) in the experimental data. 
This minimum reflects the presence of an a relaxation in the theoretical dynamics, which 
is neither observed in the rheological response (Fig. [8]) nor in the dynamic light scattering 
data (Fig. Hj). We attribute this discrepancy to the above mentioned differences between 
the dynamics of the experimental system and the glassy dynamics implicit in the theoretical 
model. 

C. Polymer Concentration Dependence of the Elastic Modulus by Mode Coupling 
Theory 

In addition to the full frequency dependence discussed above, we now consider the Cp- 
dependence of the elastic modulus G' at a fixed (low) frequency u. We compare our ex- 
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perimental data obtained for u = 0.1 rad/s (Fig. [TTl filled circles) to MCT predictions 
based on the F^'^^ model as above (open triangles), and MCT-PRISM predictions 35, 87|. 



MCT-PRISM was recently applied to silica spheres-PS mixtures at volume fraction ^ 0.4 
and shown to correctly reproduce the Cp-dependence of the measured G'{uj=l Hz) over a 
range of Cp 35|. In this study, a gel was experimentally defined on the basis of rheological 
measurements as a sample for which G'{uj=l Hz) > G" and G' > 10 Pa. Applied to our 
measurements this definition implies a gelation boundary at Cp/c* = 0.25. (In contrast, 
according to our criterion for a gel, the sample with Cp/c* = 0.25 is not considered a gel, 
because it still shows structural relaxation (Fig. HI) and thus a fluid-like response.) The range 
of polymer concentrations investigated in 35| should thus be compared to 0.25 < Cp/c* < 0.5 
in our measurements, i.e. in the vicinity of 'our' gelation boundary. For consistency with 
35| we also scale G' by 

MCT-PRISM predictions for G' (Fig. [TTl dashed line) are considerably larger than our 
experimental data. In order to compare the functional dependence we rescaled the theo- 
retical predictions. The best match with the experimentally observed trend was obtained 
for a scaling factor of about 1/4 (dotted line). This results in fair agreement; the increase 
of the experimentally observed G' is slightly more pronounced than predicted. A scaling 
factor was already proposed earlier 35| to account for the difference in the structure as- 
sumed in MCT-PRISM, a homogeneous fluid, and the silica-PS gels, which show structural 
heterogeneities in small-angle x-ray scattering experiments 65|. Structural heterogeneities 
are also present in our samples with Cp/c* > 0.25 according to our static light scattering and 
Die microscopy experiments (Figs. [H [3]). The scaling factor giving quantitative agreement 
between experimental data and MCT-PRISM predictions was calculated from the ratio be- 
tween the particle density and the density of particle clusters (with a size corresponding 
to the characteristic length of structural heterogeneities) 
of structural heterogeneities, derived from scattering experiments 
with increasing Cp for silica-PS gels. In our case, however, the characteristic length ( shows 
a strong Cp-dependence within the range of interest 0.25 < Cp/c* < 0.5 (Fig. [3]). These 
different trends of ( lead to different functional dependencies of G'. Nevertheless, rescal- 
ing the MCT-PRISM predictions by our C(cp) makes the agreement between MCT-PRISM 
predictions and our experimental data worse. 

There is some ambiguity in the determination of the gelation boundary. If we choose 



351]. The characteristic length 



651 ]. remained constant 
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FIG. 11: Elastic modulus G as a function of polymer concentration Cp/c*y G' has been normalized 
by the polymer-colloid size ratio ^, the particle volume {2R)^ and the thermal energy kB^. The 
experimentally determined G' has been extracted from dynamic frequency sweeps at the lowest 
measured frequency, u) = 0.1 rad/s, (•) and is compared to F^2^ (A) and MCT-PRISM predictions 
for G'{uj^O). The dotted hne is the MCT-PRISM prediction obtained for 4""' = 0.25c; and a 
scaling factor about 1/4. The solid line is is the MCT-PRISM prediction obtained for cf ' = 0.32c; 
and no scaling factor. The dashed line indicates the slope G' ~ (cp/c*)'^'^ 



Cp/c* = 0.32 instead of 0.25, this would possibly agree better with the definition of the 
gelation boundary by the theory, namely the transition from zero to a finite value of G. 
Moreover, this shift might be justified by differences in the structure factor between the 
silica-PS mixtures and our samples. With a gelation boundary at Cp/c* = 0.32, no scaling 
factor is needed (Fig. [IH solid line). For 0.32 < Cp/c*p < 0.5 the MCT-PRISM predictions 
describe the data well, while there are considerable discrepancies for Cp/c* > 0.5 and no 
data for larger Cp/c* are available for comparison 35|]. They might be due to the difficulty 
to treat non-equlibrium states well inside the gel region or to account for the changes in 
attraction range and strength beyond the overlap concentration c*. 

These comparisons indicate that MCT and MCT-PRISM do not entirely capture the 
mechanisms responsible for the evolution of the shear moduli approaching and entering 
the gel region. Structural heterogeneities migh cause the observed discrepancies in both 
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cases: First, structural heterogeneities can induce dynamical heterogeneities, which are not 
included in the glassy dynamics of the Fi2^ model. Second, the MCT-PRISM model is based 
on structure factors of fluid-like equilibrium structures, not heterogeneous structures. 



D. Polymer Concentration Dependence of the Elastic Modulus Inside the Gel 
Region 



The structure of a gel can be considered as closely packed fractal clusters [10|, |88|]. We 
estimate the effect of changes in this microscopic structure on the Cp-dependence of the 
elastic properties. For large clusters, the elastic behavior of the gel will be dominated by 
the deformation of clusters (strong-link regime) and the elastic constant of a cluster, K/^, is 
expected to depend on the size of its backbone, i.e. decreases with increasing cluster size. 
For small clusters, intercluster links will deform before clusters deform (weak-link regime). 
The number of particle-particle links between clusters is smaller than the average number of 
particle bonds inside a cluster, as evidenced by confocal microscopy measurements of similar 



colloidal gels 
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71| . Hence the elastic constant of the system will be dominated by the 



elastic constant o: 



G ~ K/C [88 
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intercluster links, Ki. In both cases, the total elastic contant scales as 



90|. 



Due to the large colloid volume fraction ^ 0.4, our gels consist of small clusters (Fig. [3]) 
and hence are expected to be in the weak-link regime. To obtain a scaling relation between 
G' and Cp, we have to determine the Cp-dependence of Ki and (. The correlation length ( 
sharply decreases just above the gelation boundary but then remains approximately constant 
well inside the gel region (Fig. [3]). In a minimal model, the elastic constant of intercluster 
links, Ki, is expected to depend on the number m of particle-particle contacts between 
clusters and the interaction between two particles at contact, Uq, i.e. Ki ~ mf/g. We 
assume that m does not depend on Cp. The Cp-dependence of Uq (Fig. [TOl inset) can be 
fitted by a power-law dependence, Uq ~ (cp/c*)°'^, inside the gel region. This results in 
G' ~ Ki/C ~ ttT'Uq/C Uq (cp/c*)°'^. This scaling is in agreement with G'{cp) observed 
in experiments (Fig. [TT], straight solid line). Therefore, this simple model seems to capture 
the essential mechanism leading to the elasticity of the gels, namely the intercluster links. 



28 



V. CONCLUSIONS 



We investigated the structural, dynamical and rheological properties of colloid-polymer 
mixtures with an intermediate colloid volume fraction, (p = 0.4, as a function of increasing 
polymer concentration, corresponding to increasing interparticle attraction. These samples 
covered a broad range from liquids to gels. The macroscopic gelation boundary was deter- 
mined by tube inversion. 

The structure of the samples was investigated by static light scattering and microscopy. 
Within the liquid we observed the formation of increasingly larger structures, especially 
when approaching the gelation boundary. Increasing attraction induces the formation of 
particle clusters which, at sufficiently large attraction, interconnect to form a space-spanning 
network. The maximum cluster size and maximum structural heterogeneity is observed at 
the gelation boundary. Within the gel region, increasing attraction leads to a more uniform 
structure with a reduction in the characteristic length scale. This trend is reminiscent of 
critical behavior expected for phase separation that is arrested by gelation. 

Within the entire liquid phase, the dynamics shows an ergodic response. Upon approach- 
ing gelation, the short-time (in-bond) diffusion as well as the long-time diffusion, which 
leads to the final structural relaxation of the system, slows down. The more pronounced 
slowing down of the long-time decay indicates the approach of gelation and its increasingly 
stretched exponential form suggests that heterogeneities in the density fiuctuations are in- 
creased and the distribution of length scales is broadened, reminiscent of clustering. This 
is consistent with the increasing correlation length observed in static light scattering and 
microscopy. This suggests that structural and dynamical heterogeneities, namely clustering 
and the formation of transient networks, are precursors of gel formation. 

Rheological measurements in the linear viscoelastic regime show, with increasing polymer 
concentration, a shifting of the crossing point of the elastic and viscous moduli corresponding 
to a transition from a liquid-like to a gel-like response at a characteristic frequency (time). 
This characteristic time could be related to the 'bond' lifetime estimated by the time needed 
to escape from the interparticle attraction. Bond breaking was found to be the dominating 
process close to the gelation boundary, while at lower polymer concentrations the relaxation 
appears to be related to particle or cluster diffusion. 

The elastic modulus of samples approaching the gelation boundary have been compared 
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to MCT predictions within the F^'^'^ model 50|. The predictions reproduce the depen- 
dence of the modulus on attraction strength, i.e. polymer concentration. However, the 
frequency-dependence of the modulus shows discrepancies. This could be due to structural 
and dynamical heterogeneities, which are not included in the theory. 

We compared the same data, the dependence of the elastic modulus on polymer concen- 
tration, also to MCT-PRISM predictions 35|]. Comparison between experiment and theory 
requires a consistent definition of the gelation boundary and/or a scaling factor. This is 
necessary, because, again, heterogeneities are present, which are not considered by the the- 
ory. With an appropriate choice, agreement can be obtained in a limited range of polymer 
concentrations, i.e. attraction strengths. 

In the gel region, the shear moduli show a solid-like behavior, with G' weakly frequency- 
dependent and always larger than G", which shows a minimum in the frequency dependence. 
At fixed frequency, G' increases almost linearly with polymer concentration. We suggest that 
this is consistent with a fractal model for gel elasticity in the so-called weak-link regime (simi- 
lar to low volume fraction gels 10, llSl, ISSj ) . Within this model, the almost linear dependence 
of G' results from the increase of the energy of inter-cluster links with increasing polymer 
concentration, while, at the same time, the cluster size in the gel remains approximately 
constant. This implies that, due to the heterogeneous structure of the samples, the elasticity 
of the gels is dominated by cluster-cluster links rather than particle-particle bonds. 
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